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Abstract. We present an upper limit of galactic 
diffuse gamma-ray flux, measured with the GAMMA 
experiment at energy about 175 TeV. The results 
were obtained using selection of muon poor extensive 
air showers at mountain level (700 g/cm 2 , Mt. Ara- 
gats, Armenia) for 5 GeV muon energy threshold. 
An upper limit for the differential flux at energy 
E 1 ~ 175±2o TeV is equal to (5.8 - 7.0) • 1CT 12 
{erg ■ m 2 ■ s ■ sr)™ 1 for 95% confidence level Q. 

Keywords: Cosmic rays, gamma ray, energy spec- 
trum 

I. Introduction 

Ultra-high energy (E > 100 TeV) galactic gamma- 
radiation is an important source of information about the 
origin of cosmic rays and their propagation. According 
to conventional model of cosmic rays the expected 
flux of ultra-high energy 7-rays in the energy range of 
0.1 — 1 PeV has presumably to have hadronic origin 
from sources distributed within radii ~ 1 — 0.01 Mpc 
JT] respectively. 

In the very-high (TeV) energy region 7-ray flux 
was measured by ground-based systems and single 
Cherenkov telescopes (Whipple [2|, CANGAROO 
0, MAGIC E), H.E.S.S. 0, MILAGRO [6|). The 
measurements in the ultra-high energy region (E > 100 
TeV) are still poor and made only by extensive air 
showers (EAS) arrays at Chacaltaya Q, MSU 0, 
Tien-Shan 0, CASA-MIA flj], EAS-TOP fl2l . 
KASCADE fl_D and Grapes-III ifTUl. 

This paper is devoted to measurements of the diffuse 
gamma-rays with GAMMA EAS array J_4), |Q__ at 
Aragats mountain observatory, where correlation of 
observable shower parameters with primary energy is 
about 0.98. The primary nuclei (predominantly H) and 
7-showers were discriminated using no-muon signal 
from underground muon scintillation carpet. 



II. GAMMA EXPERIMENT 

GAMMA experiment 04), EG] is the ongoing study 
of primary energy spectra in the range of 10 14 — 10 17 
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Surface EAS Array 




Fig. 1. Layout of GAMMA experiment. 



eV using EAS array with 33 concentric disposed 
scintillator stations (each of 3 x lm 2 ) and underground 
scintillation carpet with 150 scintillators (each of lm 2 ) 
to detect the shower muon component with energy 
Efj, > 5GeV/ cos 8, where 8 is the shower zenith angle. 
Layout of GAMMA facility is presented in Fig. 1. 

The detailed description of experiment, results of 
exploration of p, He, O-like, and Fe-like primary nuclei 
energy spectra derived from parametrized EAS inverse 
problem solution are presented in lfT4l . The all -particle 
primary energy spectrum obtained from event-by-event 
analysis is published in fl5l . 

Here, the discrimination of 7-showers from primary 
nuclei induced showers is performed on the basis of 
following 6 conditions: 

1) the reconstructed shower core coordinates is 
distributed within radius of R < 15m; 

2) shower zenith angles 8 < 30°; 

3) reconstructed shower size N c h > 10 5 ; 

4) reconstructed shower age parameters (s) is distributed 
within 0.4 < s < 1.5; 

5) goodness-of-fit test for reconstructed showers 
X 2 < 2.5; 

6) no-muon signal is recorded for detected showers 
satisfying the previous 5 conditions. 

The selection criteria and 7-shower discrimination 
rule (6) above were obtained using CORSIKA shower 
simulation code lfl6ll for the NKG and EGS modes 
in the frameworks of the SIBYLL 11711 interaction 
model. Simulations were done for 4 nuclear species: 
p, He, O, Fe using united energy spectral index 
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Fig. 2. Trigger efficiency of GAMMA EAS array for different primary 
particles. Dashed lines correspond to the exponential approximations: 
W trg (A, N ch ) = 1 - ot A oxp (~N ch /N (hA ), where ot A and N 0tA 
parameters depend on primary particle (A). 



Fig. 3. Expected reconstruction error (upper panel) and average bias 
(lower panel) of shower size (N c f l ) for different primary particles 
(symbols). A c ;j = ln(Af*. /AT C ^) and N*. is an estimation of shower 
size N ch . 



7 = —2.7 114L Simulated samples were equal 
7.5 x 10 5 , 10 5 , 7 x 10 4 and 5 x 10 4 for p, He, O, Fe 
nuclei and NKG mode of CORSIKA. The samples for 
the EGS mode of CORSIKA were equal to 2.5 x 10 4 
for primary 7-quanta, 7.5 x 10 4 for primary protons 
and 3 x 5000 for He,0 and Fe primary nuclei. The 
simulation strategy and reconstruction method for 
shower size (N c h), age parameter (s), core coordinates 
(xo,yo) and shower zenith angle {8) were the same as 
in HI. 

The shower trigger efficiency and shower size 
reconstruction errors (Ajv ch and <JN ch ) are presented 
in Figs. 2,3 respectively. The observed differences of 
reconstructed shower size biases Ajv ch for different 
primary particles (Fig. 3, lower panel) stems from 
differences of corresponding lateral distribution 
functions of shower particles. 

The distribution of detected shower age parameters 
(GAMMA data) in comparison with expected 
distributions for primary p, He, O, Fe nuclei are 
presented in Fig. 4 (left panel). The primary elemental 
composition and energy spectra were taken from 
solution of parametrized EAS inverse problem lfl4l 
in the frameworks of the SIBYLL interaction model 
(Fig. 5, shaded area |[T4l ) and were extrapolated up 
to the 100 TeV energy region. The reliability of this 
extrapolations stems from data fl8l . 

The right panel of Fig. 4 shows distribution of 
shower age parameters for selected no-muon signal 
events (shades GAMMA data area) in comparison 
with corresponding expected distributions from the 
simulated 7-showers and background proton showers. It 
is well seen, that the EAS age parameter is also a data 
carrier about primary particle (7-showers are younger). 
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Fig. 4. Shower age parameter (s) distribution for all showers 
(left panel) and no-muon detected showers (N^ = 0, right panel). 
Simulated data for the primary mixed composition (All nuclei, left 
panel) and primary gamma ray (right panel) are normalized to the 
corresponding GAMMA experimental data. 



However, we did not include yet the age parameter 
in the 7-showers selection criteria. Results in Fig. 4 
we use only as indication of an agreement between 
simulated and corresponding detected distributions. 

The detected muon number spectra in the 
normalization of probability density function for 
different shower size thresholds (N c h > 10 5 ,2 x 10 5 
and 4 x 10 5 ) are presented in Fig. 6 (hollow symbols) 
in comparison with the corresponding expected spectra 
from different primary particles (^f,p,He,0,Fe) and 



PROCEEDINGS OF THE 31 st ICRC, LODZ 2009 



3 




10 10 10 

g Primary particle energy, E (TeV) 

Fig. 5. Primary energy spectra and all-particle energy spectrum taken 
from fl4j (shaded area) and corresponding extrapolations to the 100 
TeV energy region. 
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Fig. 6. Normalized detected muon number (iV M ) spectra for different 
shower size thresholds (10 5 ,2 X f0 5 ,4 X 10 5 ). Hollow symbols 
(circle, square and triangle) are GAMMA experimental data. The 
symbols in 7egs. (p — F&)egs an d (p — Fz)nkg columns 
correspond to simulated data for the primary 7 and mixed composition 
(p, He, O, Fe |14| ) computed using the EGS and NKG modes of 
CORSIKA. 



different simulation modes (NKG, EGS) of CORSIKA. 
Energy spectra and elemental composition of primary 
nuclei (Fig. 5) used in the Figs. 4,6 were taken from 
|[T4l and applied for the energy region E > 100 TeV. 



A. Energy estimation 

The energy of primary particle is estimated using 
event-by-event method lfT31 according to the empirical 
expression: In E A = A\ In N c h + A2/ cos 9 + A3, where 
E is in GeV, parameters A\, A2 and A3 are determined 
using goodness-of-fit test for simulated database and 
depend on primary particle A. 



Fig. 7. Primary energy (E) and corresponding shower size (N^) 
distributions at observation level for 5000 primary protons (bold dot 
symbols) and 5000 7 (gray dot symbols). Hollow circles correspond 
to the proton showers with no-muon signal from underground muon 
carpet (jV^ = 0). Solid and dashed lines are the log-linear approxima- 
tions (see text) for primary 7 and p correspondingly. Inset histograms 
are E Pi -y/N c h distributions for primary proton and 7 (shaded area). 



The corresponding accuracies providing for x 2 — 1 
is described by the log-linear functions 
a{E A ) = E A - S A \-a(E A /10 5 ), where e = 0.22, 
0.30, 0.33 and 8 = 0.01, 0.02, 0.05 for primary 7, 
p and proton induced no-muon detected shower (po) 
respectively. 

The primary energy reconstruction efficiencies are 
presented in Fig. 7 for different primary particles. The 
lines represent the log-linear empiric expression above 
for primary proton showers (dashed line) and 7-showers 
(solid line) respectively. It is seen, that the proton pro- 
duced no-muon signal events (hollow circles) practically 
similar to 7-showers. The inset histogram shows the 
distribution of E pn /N c h ratio for p and 7 primary 
particles respectively. 



III. Gamma ray flux 

98000 shower events were selected for operation 
time T = 3970 h of GAMMA experiment. Number 
of detected shower events versus number of detected 
muons (NfJ for different shower size thresholds N c h > 
10 5 ,2 x 10 5 and 4 x 10 5 are presented in Fig. 8 
(histograms with shaded area). The symbols in Fig. 8 are 
the corresponding expected number of events simulated 
using the CORSIKA code for primary energy spectra 
lfl4l presented in Fig. 5. The simulations were carried 
out for two modes of CORSIKA to get high accuracy 
of simulation (EGS mode) and large simulated sample 
(NKG mode). 

The agreement of simulated and detected muon spec- 
tra in all measurement region and lack of statistically 
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Fig. 8. Detected (histogram lines) and expected (symbols) muon 
number (N^) spectra for different shower size thresholds (10 5 ,2 X 
10 5 ,4 X 10 5 ) and different mode (NKG, EGS) of CORSIKA. 



significant excess of no-muon signal events (Fig. 8) 
allowed us to estimate only an upper limit of 7-ray flux 
according to the expression 



J = 



2y/M /{W llNlt =oWtr gn ) 1 

SnTcosO ' AE 



(1) 



where Mo is the number of no-muon detected showers, 
W 7 ,iv„=o lS m e probability to detect no-muon signal 
for 7-showers (see Fig. 6), Wtrg(E~) is the trigger 
efficiency (see Fig. 2), cos 8 = 0.94 is the average 
shower zenith angle, S and are the EAS core 
detection area and corresponding solid angle. 

The obtained upper limit of differential 7-ray flux in 
the energy range 100-300 PeV are presented in Fig. 9 
(black downward triangle symbol) in comparison with 
CASA-MIA flU, KASCADE 03 and EAS-TOP ifHl 
data. 



IV. Conclusion 

An upper limit of 7-ray differential flux at energy 
E 7 ~ 175+20 TeV obtained with GAMMA experiment 
is equal to (5.8 — 7.0) • 10~ 12 (erg ■ m? ■ s ■ sr)^ 1 for 
95% confidence level and it is in close agreement with 
the CASA-MIA data ifTTIl . 

The lower limit for the primary energy spectra and 
elemental composition obtained with the GAMMA 
experiment lfT4ll can be extended to the lower energy 
region up to about 100 TeV energies. 

We are going to increase the underground muon 
carpet area up to 250 m 2 to improve 7/proton showers 
discrimination efficiency. 



Fig. 9. Upper limit of gamma ray flux derived from detected no- 
muon showers (black downward triangle symbol). The gray symbols 
are the CASA-MIA [TT], KASCADE QH and EAS-TOP Qz) data 
taken from 1131 . The lines are expected Galactic diffuse background 
flux from CD- 
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